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ABSTRACT 

Water is a crucial molecule in molecular astrophysics as it controls much of the 
gas/grain chemistry, including the formation and evolution of more complex or- 
ganic molecules in ices. Pre-stellar cores provide the original reservoir of material 
from which future planetary systems are built, but few observational constraints 
exist on the formation of water and its partitioning between gas and ice in the 
densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, 
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we report on the first detection of water vapor at high spectral resolution toward 
a dense cloud on the verge of star formation, the pre-stellar core L1544. The 
line shows an inverse P-Cygni profile, characteristic of gravitational contraction. 
To reproduce the observations, water vapor has to be present in the cold and 
dense central few thousand AU of L1544, where species heavier than Helium are 
expected to freeze-out onto dust grains, and the ortho:para H 2 ratio has to be 
around 1:1 or larger. The observed amount of water vapor within the core (about 
1.5 xl0~ 6 M Q ) can be maintained by Far-UV photons locally produced by the im- 
pact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the 
icy mantles, liberating water wapor in the core center. Our Herschel data, com- 
bined with radiative transfer and chemical/dynamical models, shed light on the 
interplay between gas and solids in dense interstellar clouds and provide the first 
measurement of the water vapor abundance profile across the parent cloud of a 
future solar-type star and its potential planetary system. 

Subject headings: Astro chemistry — Line: profiles — Stars: formation — ISM: 
clouds — ISM: molecules 



Introduction 



The ubiquity and importance of w ater in the Universe has become clear after the launch 
of the Herschel Space Ob servatorjj] ( Pilbratt et al. 2010 ). Water vapor has been found 



ies and active galactic nuclei ( Ivan der Werf et al. 



aroun d moons and comets (lHartogh et al.ll2011a 
2011 ). star forming regions Iceccarelli et al. 



bh, in protoplanetary disks (jHogerheijde et al. 



2010: Ivan Dishoeck et al. 



2010 



Wcifi et al. 



20111 ). starburst galax- 



20101 ). One important 



class of objects is missing from this list: the centrally concentrated cores of molecular clouds, 
which precede the formation of stars and planetar y systems, as previou s studies have onl y 
provided upper limits of their water vapor conte nt (Bergin fc Snellll2002l : iKlotz et al.ll2008l ). 
These so-called pre-stellar cores, as defined by Crapsi et al. ( 20051 ). are centrally concen- 
trated, with central densities > 10 5 cm -3 , they have central visual extinctions >50mag, 
large degrees of CO freeze-out (> 90%) and deuterium fraction (>20%), and they show clear 
signs of contraction motions in high-density tracer spectra. Pre-stellar cores are typically 



cold (T < 20 K) an d detailed studies of two pre-ste 



of only 7K: L1544 (ICrapsi et all 120071 ) and L183 (IPagani et al.l 120071 ) 



lar cores have foun d central temperatures 



1 Herschel is an European Space Agency space observatory with science instruments provided by European- 
led principal investigator consortia and with important participation from NASA. 
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Pre-stellar cores represent the initial conditions in the process of star and planet forma- 
tion. Thus, a direct measurement of water toward a pre-stellar core would provide a unique 
opportunity to unveil the amount and the chemistry of water at the dawn of a new solar-type 
system. Water has been clearly detected in solid form in molecular clouds, where pre-stellar 



cores are embedded, or in dense cores hosting young stellar objects (IWhittet fc Duleylll991 



Boogert et al.ll2008l ; IChiar et al.ll201ll ). However, these near- and mid- infrared solid-state ob- 
servations cannot be extended to the realm of pre-stellar cores, given their large extinctions 
and starless nature. How can we measure the water abundance in a pre-stellar core? 

The cold dust present in pre-stellar cores has its peak of emission at wavelengths > 
290 /iin and dust continuum emission is expected to be detectable at 557 GHz, the frequency 
of the ground state transition of ortho-H 2 (Jr a — Jk c = lio-loi)- Thus, the 557 GHz water 
line (difficult to be collisionally excited in cold gas) could be se en in absorpti o n aga inst the 
far-infrared background produced by the core central regions. ICaselli et al.l (120101 ) indeed 
detected an absorption feature at a 5 a level using the wide-band spectromete r (WBS) of 
the Heterodyne Instrument for the Far Infrared (HIFI; Ide Graauw et al.l 120101 ) toward the 
pre-stellar core L1544. However, the spectrum did not have enough sensitivity and spectral 
resolution to study the line profile and allow quantitative conclusions. Here we report the re- 
sults of a deeper search with 9 times greater spectral resolution and 2 times better sensitivity 
that has provided, for the first time, a velocity-resolved spectrum of water towards L1544, 
revealing the presence of water vapor in the central few thousand AU (the core center) for 
the first time. The technical details of the observations are described in Section [21 the data 
are presented in Section [21 a closer look to the water line profile is given in Section 0] and 
the data interpretation is described in Section [51 



2. Observations and data processing 



As part of the WISH (Water In Star forming regions with Herschel) k ey project (Ivan Dishoeck et al. 



20111), the ortho-H 2 0(li -loi) line (frequency 556936.002±0.050 MHz; IPickett et al. 
was observed with Herschel towards L1544 (coordinates «2ooo — 05 /l 04 m 17 s .21, 5~ 



200 



2000 



25°10'42".8) for 3.3 hours on March 20, 2010 and for 10.3 hours on April 2, 2011. The data 
presented here are archived at the Herschel Science Archival, under OBSID 1342192364 and 
1342217688. The half-power beam width (HPBW) of the telescope at this frequency is 40". 
The technical setup in the two periods included dual beam switch mode, with the continuum 
optimization and fast chop options, and simultaneous use of the wide-band (WBS) and the 



2 http://archives.esac.esa.int/hda/ui/ 
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high-resolution (HRS) spectrometers in both H and V polarization. The individual spectra 
were redu ced using the Herschel Interactive Processing Environment (HIPE, version 6.1.0; 



Ott||2010j). Subsequent analysis of the data was performed with the Continuum and Line- 



Analysis Single dish Software (CLASS) within the GILDAS packag^f): the line intensity units 
were converted from an tenna to brightness temperature (7b) using a beam efficiency of 0.76 
(IRoelfsema et al.l 120121 ); the H and V spectra were analyzed individually and found similar 
within the noise (therefore, the two spectra were averaged). After binning to the nominal 
spectral resolution (0.59 km s -1 for WBS and 67ms -1 for HRS), the root-mean-square (rms) 
noise level is 0.93 mK in the WBS spectrum and 3.8 mK in the HRS spectrum. The dust 
continuum emission flux at 557 GHz (10.2±0.2mK) was obtained by dividing by 2 the in- 
tensity of the WBS 2011 dual beam switch spectrum, degrading its spectral resolution by a 
factor of 50 and measuring the mean and standard deviation of the brightness temperature 
of the remaining channels. 



3. The water spectrum toward L1544 

The (continuum-subtracted) spectrum of ortho-H20(lio-loi) observed with Herschel/HIFI 
toward L1544 is shown in Fig. 1. In the t op figure, the low spectral resolution HIFI- 



WBS spectrum of ortho-H 2 taken in 2010 (jCaselli et al.ll2010l ) is superposed on the final 



(2010+2011) spectrum, confirming the absorption feature. The integrated intensity of the 
feature is now 10 times the associated error. The HRS spectrum smoothed at the WBS ve- 
locity resolution is also shown. The WBS absorption feature has a width at half-maximum 
of 2kms _1 , an optical depth of 1.0 ±0.2 and a column density of ortho-H 2 molecules in the 
ground state of (2.0+0.4) xlO 13 cm -2 . This is the first accurate water vapor column density 
estimate in a dark cloud, away from active sites of star formation. However, this water does 
not exclusively belong to L1544. In fact, the WBS absorption feature spans a large range of 
velocities and the bulk of the absorption appears shifted by about +2 km s" 1 compared to 
the local standard of rest (LSR) velocity of L1544, Vlsr = 7.2 km s -1 , marked by the dotted 
vertical line. To unveil the water content in L1544, velocity-resolved analysis is needed. 

The yellow strip in Fig. la encloses the velocity range displayed in Fig. lb, where the 
final (2010+2011; red) and 2010 (black) HRS spe ctra of ortho-H 2 Q ( Iiq-Iqi) are superposed 



on the CO (J = 1— > 0) spectrum (grey histogram; ICaselli et al.ll2010l ). The intensities of the 
2010 water line and the CO line have been reduced by a factor of 2 and 100, respectively, 
for displaying purposes. The final HRS spectrum reveals three features not noticeable in 



3 http://www.iram.fr/IRAMFR/GILDAS. 
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the upper panel: (1) an emission feature shifted by -Oikrns" 1 compared to the Vlsr of 
L1544 (dotted vertical line); (2) an absorption feature centered at the Vlsr of L1544; (3) 
an absorption feature +1.9 km s -1 away from the Vlsr of L1544 and close to the secondary 
feature seen in CO (1-0), which is probably an independent cloud along the same line of 
sight (not to be considered further in this study as physical characteristics of this cloud are 
not yet known). The results of independent Gaussian fits to the three spectral components 
are listed in Tab. 1. These features are spectrally diluted in the WBS spectrum, with the 
emission+absorption (components 1+2 in Fig. lb) toward L1544 resulting in the shallower 
absorption feature seen in the WBS spectrum at the velocity of L1544. The HRS spectrum 
thus provides information about the distribution of water along the line of sight. 

The water absorption feature seen at the velocity of L1544 can be used to estimate the 
amount of water vapor present within the pre-stellar core. However, the intensity of the 
absorption is consistent with saturation (Tab. 1), so only a lower limit can be gauged. From 
a curve-of-growth analysis, the optical d epth of the LI 544 absorption feat ure is r > 1.5. 
Assuming a H2O ortho-to-para ratio of 3 (ILis et al.ll2010i ; iHerpin et al.ll2012l ) and a width at 
half-maximum of 0.5 km s -1 (Tab. 1), the total (ortho+para) H 2 column density across the 
pre-stellar core is >1.0xl0 13 cm -2 , corresponding to a fractional abundance (wi th respect to 



H 2 m olecules) >1.4xl0 10 (using an H 2 column density iV(H2) — 7xl0 22 cm 2 ; ICaselli et al. 



20101 ). A detailed radiative transfer analysis is needed to unveil the amount of water within 
L1544 (see Section EJ). 



Central water inflow 



To better understand the nature of the water emission in L1544, the water spec- 
trum between_5^and_8. 5km s _1 is plotted in Fig. 2 together with the o rtho-H 2 D + (l-in-lii ) 
line ( ICaselli et al.ll2003l ). This line ha s a critical density of 10 5 cm -3 (jHugo et al.l I20091 ) 



traces gas at radii of ab out 5000 AU (jVastel et al.l I20061 ) and it is only moderately thick 
(Ivan der Tak et al.ll2005l ). The ortho-H 2 0(lio-loi) line has a significantly larger critical den- 
sity (>1.7xl0 7 cm -3 , for an ortho-to-para H 2 ratio < 3), and the line is highly subthermally 
excited in the region enclosed by the Herschel beam . These are ideal conditions to observe 
the characteristic "infall asymmetry" (jEvansI Il999l ): foreground cloud material, while col- 
lapsing, is absorbing the high-velocity portion of the line profile, producing a brighter blue 
peak. In the presence of absorption against the continuum background emission, this is 
equivalent to an inverse P-Cygni profile, as observed in L1544. The water blue emission fea- 
ture reveals the inner backside of the core, which is approaching the cloud cen ter. Previous 



signatures of infall of pre-stellar cores were found using the N 2 H + (l-0) line ( ICrapsi et al. 
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20051 ). which traces gas with densities ~10 5 cm -3 . The water line is tracing the gravitational 
contraction of gas at densities larger than 10 6 cm -3 , which resides within a radius of 1000 AU 
(see Sect. \5jj. Inverse P -Cygni profiles have also been found in 6 out of 15 Class objects 
( IKristensen et al.ll2012l ). highlighting water's ability to pick up infall motion along the line 
of sight. 

The detection of emission (component #1 in Fig. lb) is crucial for our understanding of 
where water resides in this system. From chemical considerations alone it has been difficult to 
posit the presence of water in the dense center. We know that water ice is prese nt in the less 
dense regions of the extended Taurus Molecular Cloud ( IMurakawa et al.ll2000l ). accounting 
for almost 60% of the total available gas phase Oxygen (3.2xl0~ 4 with respect to H nuclei; 
Meyer et al.lll998l ). Such a large amount of water cannot be produced via gas-phase routes 
and then deposited on top of dust grai ns, as gas-phase chemistry can only account for about 
one thousandth of the ice abundance ( jLee et al.lll996l ). Thus, it is now well-established that 
water ic e is formed on the su rface of dust particles via hydrogenation of atomic and molecular 
oxygen ( lloppolo et al.l 120081 ) . 



It is also well-established that atomic and molecular species start to freeze-out onto dust 
grain s in pre-stellar core s at volume densities larger than a few times 10 4 cm -3 ( IWalmsley 
1 99 ll : ICaselli et al.lll999l ). L1544 is among the densest of these objects, with central values 
>10 6 cm~ 3 (Sect. [5]). The relatively volatile CO molecules (evaporation temperature of 
~20 K) are mostly frozen onto grains in the core centers. Water, with an evaporation 
temperature of ~ 100 K, should be more difficult to remove from dust grains and is therefore 
less likely to be observed via emission in the gas-phase. However, the detection of emission 
in our spectrum tells a different tale. A small amount of water vapor must be present in the 
core center, contrary to our theoretical expectations, along with its presence in the outer 
layers seen in absorption (component #2 in Fig. lb). 



5. The importance of cosmic-rays 



The chemistry, dynamical motions, and excitation are intimately intertwined in this 
object, resulting in saturated absorption and emission. This requires detailed physical, 
chemical and radiative transfer models to more accurately estimate the water content of 
L1544 and elucidate the key aspects of the velo city field in this object. Current state-of- 
the-art chemical models ( IHollenbach et al.ll2009l ) posit that small amounts of water vapor 
exist (containing less than a percent of the available cosmic oxygen) in dense regions that are 
well below the evaporation temperature, when exposed to the far-ultraviolet (FUV) starlight 
pervading interstellar space. When present, UV photons photodesorb water molecules and 
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other species from the icy mantles (jOberg et al.ll2009l ). Thus, the desorption depends on 
the intensity of the interstellar FUV field, which drops to tiny levels within dense clouds. 
There is a small internal radiation field in these regions induced by cosmic rays interacting 
with H 2 , leading to a dim FUV field about 1000 times less intense than th e interstellar field 
( 
( 



radiative transfer code MOLLIE flKetolll990l : Keto et all 12004 ; iKeto fe Rvbickill2010l ). 



Gredel et al.lll989l) . We have placed these mechanisms withi n the model of IHollenbach et al. 
20091 ). which we follow here, in ta ndem with our dynamical (Keto fc Caselrill2008 . 2010 ) and 



Besides the comprehensive chemical network, we also consider a simplified oxygen chem- 
istry to be easily coupled with our hydrodynamical code. The details of this chemical code 
will be presented in a separate paper (Keto et al., in preparation). Here we briefly mention 
that the formation of water is assumed to happen on dust grain surfaces by the reaction of 
OH and H. Gas phase routes to the formation of water in cold gas are not included. We allow 
gas-phase O and OH, as well as H 2 0, to freeze-out onto dust grains, but assume that the 
frozen O and OH are immediately converted into H 2 ice. CO is treated independently from 
the oxygen chemistry. Photodesorption (including that due to the cosmic-ray induced FUV 
field), photodissociation and desorption due to dire ct cosmic-ray imp a cts w ith dust grains 
are also taken into account in the model, following IHollenbach et al. J2009h prescr i ptions . 
We adopt a cosmic-ray ionization rate of lx 10~ 17 s -1 , as deduced by lKeto fc Casellil (120 lOl ). 



and an average gr ain-surface cross section <th=1-4x 10 cm , a factor of 1.4 lower than the 
value adopted by IHollenbach et al.l (120091 ) . Overall, the difference between the simple and 
comprehensive model is within a factor of 2. 

Fig. 3 (bottom panel) shows the physical structure of L1544, as deduced from a de- 
tailed comparison between observational data (N 2 H + (l-0), C 18 O(l-0), C 17 O(l-0) line profiles 
and the temperature structure deduced from NH 3 ) and simulated observations of a dynami- 
cal/chemical model_of_asph^^ near hydrostatic equilibrium undergoing gravitational 
contraction (IKeto fc Casellil |2010| ). The middle panel displays three water vapor fractional 
abundance (with respect to H 2 ) profiles (color curves): (1) the red profile is the result of 
the full chemical network applied to the physical structure of L1544, without the internal 
FUV field produced by cosmic rays; (2) the green curve shows the water profile obtained in 
a self-consistent chemical-dynamical model using the simplified chemistry, without the inter- 
nal FUV field (note that results are similar to the more comprehensive model); (3) the blue 
thick curve is the result of the simplified chemistry with the inclusion of the dim FUV field 
produced by cosmic rays. The black curve is the water ice abundance (divided by 10000) 
calculated with the full chemical network inclusive of the internal FUV field. The first thing 
to note is that the FUV produced by cosmic-ray impacts is the main desorption agent in 
the core center, water molecules are partially photodesorbed from the icy dust mantles and 
the water vapor abundance jumps by 3 orders of mag nitude to > 10" 9 within 5000 AU. At 
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larger radii, the visual extinction (Ay) drops below 5 magnitudes, allowing interstellar FUV 
photons to photodesorb icy mantles and to maintain the water abundance around 10~ 7 . 
Beyond 20000 AU, Ay becomes less than 2 mag and the interstellar FUV photons efficiently 
photodissociate water, as well as other molecular species, explaining the drop in the water 
abundance. 

The simulated water lines obtained with the MOLLIE radiative transfer code, applied 
to the model cloud and water abundance profiles just described, are in the top panel of 
Fig. 3. The MOLLIE code has been adapted t o treat extremely sub-thermally excite d lines 
(such as the observed ortho-H^O line; see also lSnell et al.ll2000l ; iPoelman et al.l 120071 ). The 
blue, green and red line profiles refer to the corresponding H2O abundance profiles in the 
central panel. The brighter blue (dashed) profile is obtained if a 3:1 ortho:para ratio of H 2 
molecules is assumed as collision partners of ortho-H 2 0, while the fainter blue emission profile 



is obtained with a 1:1 or tho:para H2 ratio (using collisional coefficients from iDubernet et al. 



20091 ; iDaniel et al.ll201ll) . Lower ortho:para ratios (expected in regions w ith large amount of 
freeze-out; iPagani et al.l 120091 ; iTroscompt et al.l 12009c iParise et al.l 1201 it ) are not consistent 
with the present observations, unless a higher abundance of water vapor is present in the gas 
phase at high densities. Considering the simplicity of the physical model, the agreement with 
the data is remarkable. The similarity between the thick-blue curves and the data, including 
the continuum level, confirms that the observed water line prole is indeed the result of water 
molecules emitting in the core center and absorbing in the envelope. The observed extra- 
emission at low velocities suggests the presence of higher velocity motions, probably due to 
larger infall velocities within the central 1000 AU than predicted by our models (Fig. E]). 
This is the first tracer of dense gas known to be present near the core center that shows signs 
of gravitational contraction. 

In summary, Herschel has revealed the inflow of water vapor in a pre-stellar core and the 
key phenomena at work: (1) fractional abundance of water vapor >10~ 9 toward the central 
few thousand AU (the core center). This is possible thanks to the cosmic-ray induced FUV 
field, which illuminates the icy dust mantles, allowing H2O molecules to evaporate. (2) 
Volume densities >10 6 cm -3 in the core center. This allows (i) to collisionally excite the 
observed ground state rotational transition of ortho-H 2 and produce observable emission, 
and (ii) to provide large enough dust continuum emission in the FIR to allow detection of 
absorption features. (3) Gravitational contraction motions, which broaden the water line 
and allow the blueshifted emission to avoid the foreground "water absorption screen" . From 
the water vapor and ice abundance profiles deduced from our chemical models (FigJH]), the 
total amount of water vapor within the central 10000 AU is ~1.5x 10 -6 M Q , while a reservoir 
of about 2.5xlO~ 3 M (2.6 Jupiter masses) is available in the form of water ice. 
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HIFI has been designed and built by a consortium of institutes and university de- 
partments across Europe, the United States and Canada under the leadership of SRON 
Netherlands Institute for Space Research. 

Facilities: Herschel Space Observatory 
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Fig. 1. — (a) The Herschel HIFI-WBS spectrum (velocity resolu tion = 0.59 km s -1 ) of ortho- 
H20(lio-loi) observed toward L1544 in 2010 (black histogram; ICaselli et al.l 120101 ) and the 
final spectrum which combines 2010 and 2011 data (red histogram). The HRS spectrum 
smoothed at the velocity resolution of WBS is also shown (grey histogram). The dotted 
line marks the Vlsr of L1544 (7.2kms _1 ). (b) HIFI-HRS spectrum (velocity resolution = 
67ms- 1 ) of ortho-H 2 0(li -loi) (red: 2010+2011 data; black: 2010 data, intensity divided 
by 2). The grey histo gram is the CO (1- 0) line (intensity divided by 100), observed along 
the same line of sight ( ICaselli et al.ll2010f ). Gaussian fit results and integrated intensities of 
the three labelled features are listed in Table 1. The continuum has been subtracted from 
all the spectra. 
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Fig. 2. — HIFI-HRS spectrum of ortho- HgOflm-lm) (blue h istogram) compared with ortho- 
H 2 D + (l 10 -ln) (intensity divided by 30; ICaselli et al.ll2003l ) toward L1544. The water line 
shows an inverse P-Cygni profile, with emission only at blueshifted velocities, whereas the 
water in the bulk of the core is absorbing the far-infrared dust continuum emission produced 
by the core central regions. 
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Fig. 3. — Bottom panel: Physical structure of L1544 from iKeto &: Casellil ( 120101 ): gas tem- 
perature, H 2 number density and velocity profiles. Central panel: fractional abundance of 
water vapor (w.r.t. H 2 , color curves) as a function of core radius. Red curve: comprehensive 
model without the cosmic-ray induced FUV field. Green curve: simplified chemistry with- 
out the cosmic-ray FUV field. Blue curve: simplified model with the cosmic-ray FUV field. 
Black curve: fractional abundance of water ice (divided by 10000) calculated with the full 
chemical code, inclusive of the cosmic-ray FUV field. Top panel: ortho-H 2 0(li -loi) HRS 
continuum-subtracted spectrum (the horizontal black line displays the subtracted continuum 
level), compared with the results from the radiative transfer modeling, using as input the 
corresponding H 2 abundance profile in the central panel. The two blue curves refer to two 
different values of the H 2 ortho:para ratio, 3:1 (dashed) and 1:1 (thick). 
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Table 1: Gaussian fit results for the three features of ortho-H 2 0(li -loi)- 

Component Area V LS r Av T b 

(mKkms -1 ) (kms -1 ) (kms -1 ) (K) 

1 +4.7±0.8 6.79±0.02 0.24±0.05 +19±4 
+5.0±0.7 a ' 1 

2 -6.1±1.3 7.18±0.04 0.45±0.13 -13±4 
-4.7±0.7 a ' 2 

3 -9.4±1.7 9.05±0.10 1.12± 0.25 -8±4 
-6.8±l.l a ' 3 

( a ) Intensity integrated between channels with |T B | > 4mK: 
W from 6.52 to 6.93 km s" 1 , 

( 2 ) from 6.98 to 7.45 km s" 1 , 

(3) from 8.69 to 9.83 km s" 1 . 



